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It has recently become possible to form molecules in ultracold gases of trapped alkali metal 
atoms. Once formed, the molecules may undergo elastic, inelastic and reactive collisions. Inelastic 
and reactive collisions are particularly important because they release kinetic energy and eject atoms 
and molecules from the trap. The theory needed to handle such collisions is presented and recent 
quantum dynamics calculations on ultracold atom-diatom collisions of spin-polarised Li + L12 , Na + 
Na2 and K + K2 are described. All these systems have potential energy surfaces on which barrierless 
atom exchange reactions can occur, and both inelastic and reactive rates are very fast (typically 
fcinci > 10 -10 cm 3 s _1 in the Wigner regime). 

PACS numbers: 34.20.-b,34.20.Mq,34.50.Ez,34.50.-s,82.20.Ej,82.20.Kh 
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INTRODUCTION 

There have been enormous recent advances in our abil- 
ity to produce and trap samples of cold molecules (below 
1 K) and ultracold molecules (below 1 mK). Molecules 
such as NH3, OH and NH have been cooled from room 



molecular Bosc-Einstein condensates have beenproduced 
for dimers of fermionic alkali metal atoms [aJS E^, 
and the first signatures of ultracold triatomic [llj and 
tetraatomic LL2I molecules have been observed. 



Cold and ultracold molecules have many possible appli- 
cations. High-resolution spectroscopy on cold molecules 
may allow the measurement of fundamental physical 
properties such as the electric dipole moment of the elec- 
tron J^3j, the energy differences between enantiomers 
[I4I llo| and the time-dependence of the fine-structure 
constant |l6| . In addition, since molecules have a much 
richer energy level structure than atoms, they offer many 
new possibilities for quantum control. Perhaps most im- 
portantly, dipolar molecules interact with one another 
much more strongly and at longer range than atoms. 
Dipolar molecules have been proposed as qubits for quan- 
tum computers and dipolar quantum gas es are pre- 
dicted to exhibit a range of novel features |18j . 



In a recent article ja| , we reviewed the current state of 
the art of molecule production in ultracold atomic gases. 
Other authors have reviewed the cooling of molecules 
from near room temperature 0, and the theory of col- 
lisions of such directly cooled molecules [ltA |2fj, |2lJ . The 
present article is complementary to these and focusses on 
recent theoretical work on the collisions of alkali metal 
dimers formed in ultracold gases. 
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EXPERIMENTS ON ULTRACOLD MOLECULE 
FORMATION AND COLLISIONS 

There are two main methods used to form molecules in 
ultracold atomic gases. In photoassociation a pair 

of atoms undergoes a spectroscopic transition from an 
unbound atomic state very close to threshold to a bound 
molecular state. In magnetoassociation, also known as 
Feshbach resonance tuning |5j,|7j, the transition is accom- 
plished by guiding the atom pair adiabatically across an 
avoided crossing between an unbound state and a molec- 
ular state. 

Feshbach resonance tuning always produces molecules 
in very highly excited states, usually the highest vibra- 
tional state that exists in the diatomic potential well. 
Photoassociation also usually produces molecules in very 
high vibrational states, because they are the ones that 
have the best Franck-Condon overlap with the free-atom 
states. 

Ultracold molecules are initially formed in the pres- 
ence of ultracold atoms and can collide with them. For 
molecules in vibrationally excited states, there is the pos- 
sibility of vibrationally inelastic collisions, 

M 2 (»+M — >M 2 {v' <v)+M, (1) 

where v is the vibrational quantum number. Since the 
trap depth is usually much less than 1 K, such collisions 
always release enough kinetic energy to eject both colli- 
sion partners from the trap. If the molecular density is 
high, there is also the possibility of inelastic molecule- 
molecule collisions, 

M 2 (v) + M 2 (v) — ► M 2 (V < v) + M 2 (v" < v). (2) 

Molecules are not destroyed in inelastic collisions, but 
they are lost from the trap and are no longer ultracold. 

The initial experiments on molecule formation by Fes- 
hbach resonance tuning worked with bosonic isotopes 
of alkali metals [H Hi . They found that the 

molecules were lost from the trap on a timescale of mil- 
liseconds. The loss was attributed to vibrationally inelas- 
tic atom-molecule collisions with relaxation rates around 
10~ 10 cm 3 s _1 . For the case of 85 Rb 2 0, recent work 
has suggested that the loss may in fact be due to sponta- 
neous molecular dissociation by collisionless spin relax- 
ation jUizI]. However, for the other systems [H 0,111] 
the molecules are formed in truly bound states and can- 
not decay without collisions. Mukaiyama et al. jll| have 
recently measured the trap loss rate for 23 Na 2 molecules 
formed by Feshbach resonance tuning and obtained an 
atom-molecule rate coefficient fci oss = 5.1 x 1CP 11 cm 3 
s _1 for molecules in the highest vibrational state. 

Fermion dimers formed by Feshbach resonance tuning 
are a very special case. In mid-2003, four groups inde- 
pendently reported within a very short time that dimers 



of fermionic 6 Li jH US EJ and 40 K H3 could be re- 
markably stable to collisions. Cubizolles et al. [3(| and 
Jochim et al. [sil ] showed that the lifetime was partic- 
ularly large close to a Feshbach resonance, where the 
scattering length is large and positive. By the end of 
2003, three different groups H,M had succeeded in 
creating long-lived molecular Bosc-Einstcin condensates 
of fermion dimers. 

Petrov et al. [33I I34I ] analysed the stability of fermion 
dimers in terms of the long-range form of the wavefunc- 
tion. In the case where the atom-atom scattering length 
a is much larger than the range of the atom-atom po- 
tential r e , they showed that both atom- molecule and 
molecule-molecule inelastic collision rates are suppressed 
by Fermi statistics. However, their derivation applies 
only to molecules that are in long-range states, with 
a wavefunction that depends on the scattering length, 
x(r) ~ exp(— r/a). As will be discussed in more detail 
below, Cvitas et al. |35l | have shown computationally that 
there is no systematic suppression of the atom-molecule 
inelastic rate for fermion dimers in low-lying vibrational 
levels, even when a is large and positive. 

Relaxation processes have also been studied for 
molecules formed by photoassociation. Wynar et al. pjjf 
formed 87 Rb2 molecules in the second-to-last vibrational 
level of the ground excited state by stimulated Raman 
adiabatic passage (STIRAP). They estimated an upper 
bound of k\ oss = 8 x 10~ n cm 3 s _1 due to inelastic atom- 
molecule collisions. Staanum et al. investigated in- 
elastic collisions of rovibrationally excited CS2 ( 3 S+) in 
collisions with Cs atoms in two different ranges of the 
vibrational quantum number v by monitoring trap loss 
of CS2. They obtained atom- molecule rate coefficients 
close to 1.0 x 10~ 10 cm 3 s _1 for both v = 4 to 6 and 
v = 32 to 47. Zahzam et al. |38| carried out similar work 
for different rovibrational states of 3 S+ , and also consid- 
ered molecules in the 1 E+ state and molecule-molecule 
collisions. They obtained rate coefficients of 2.6 x 10~ n 
cm 3 s _1 and 1.0 x 10~ n cm 3 s _1 in the atom-atom and 
atom-molecule cases respectively, both with quite large 
error bounds. 

Because of the collisional losses for excited vibrational 
levels, intense efforts are under way to produce ultracold 
alkali metal dimers in low-lying levels (and ultimately 
in the ground vibronic state). The process of transfer- 
ring molecules from atomic or near-dissociation molecu- 
lar states, with probability density at long range, to low- 
lying states, with probability density near the diatomic 
equilibrium distance r e , is sometimes called r-transfer. It 
can in principle be achieved either by using man y p ho- 
tons to accomplish the transfer in several stages |39l or 
by using tailored ultrafast laser pulses 0> EH E3- ^ n 
favourable cases, and particularly for heteronuclear al- 
kali metal dimers [43J, |44j , it may be possible to produce 
molecules in their vibrational ground state by two-photon 
processes via excited states with mixed singlet and triplet 
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character. Sage et al. 43] have recently succeeded in cre- 
ating ultracold RbCs molecules (T w 100 fiK) in their 
vibronic ground state using a 4-photon process. The 
production rate in the current experiments is only 500 
molecules/s, but work is under way to increase it. 

Inspired by the recent progress in experimental tech- 
niques for producing and studying ultracold molecules, 
we have collaborated with Jean-Michel Launay and 
coworkers at the Universit y of Rennes to produce a series 
of theoretical studies E3, E3, EH E3| on 

interactions and collisions between spin-polarized alkali- 
metal atoms and molecules. In particular, ultra-low- 
energy collisions were studied for Na + Na2 |iR lijj , Li 
+ Li 2 (isotopically homonuclear [3^, 0] and heteronu- 
clear |B3)j an d K + K 2 [EH- Collisions between alkali- 
metal dimers and atoms present several new theoretical 
challenges that are not present for collisions of stabler 
molecules. The remainder of this review will focus on de- 
scribing this work and drawing general conclusions from 
it. 



GENERAL ASPECTS OF ALKALI METAL 
DIMER COLLISIONS 

The atoms in an ultracold quantum gas of S-state 
atoms are in states labelled by the electron spin s, the 
nuclear spin i, and the total angular momentum /. Since 
there is usually a magnetic held present, the states are 
also labelled by to/ , the projection of / along the held. Of 
particular interest here are spin- stretched states, where 
/ = /max = i + s and |m/| = /. From a theoretical 
point of view, collisions of atoms and molecules in spin- 
stretched states are simpler than others, because they 
take place entirely on high-spin potential energy sur- 
faces. For alkali metal atoms with s = ^, these are triplet 
curves (S = 1) for atom-atom collisions and quartet sur- 
faces (S — |) for atom-molecule collisions and 3-body 
recombination. For atoms in non-spin-stretched states, 
singlet dimer curves and doublet trimer surfaces are also 
required. Our work so far has focussed on spin-stretched 
states. 

The lowest quartet state of the alkali metal trimers, 
which correlates with ground-state atoms ( 2 S) and 
molecules in their lowest triplet state ( 3 S+), is desig- 
nated (I 4 A'). At first sight it might be expected that 
three parallel electrons in s orbitals would not form sig- 
nificant chemical bonds. If this was the case, then pair- 
wise additivity would be a good approximation and the 
potential energy surface would be given by 

3 

V(r 12 ,r 2 3,ri 3 ) w ^ Vj 

)■ (3) 

i<j 

The potential energy curve for the lowest triplet state 
is reasonably well known for many of the alkali metal 



dimers Pairwise-additive model poten- 

tials have been extensively used in theoretical studies of 
3-body recombination [53, [53, [5JJ, [6fJ • 

As will be seen below, pairwise additivity is actually 
quite a poor approximation for the alkali metals. Never- 
theless, it gives some insights into the energetics of atom- 
molecule collisions. If the dimer well depth is e at an 
atom-atom distance r e , then for a pairwise-additive sur- 
face the energies of various important arrangements of 3 
atoms are 

• atom M well-separated from diatom M2 with bond 
length r e , 

V(r) = V(r e ,oo, 00) = — e 

• linear trimer M3 with bond length r e , 

y(r) = T/(r e ,r e ,2r e )«-2e 

• triangular trimer M3 with bond length r e , 

V(r) = V(r e ,r e ,r e ) = -3e. 

The topology of the lowest quartet potential energy sur- 
face of an alkali-metal trimer is thus quite simple; the 
global minimum is at an equilateral triangular config- 
uration (point group D^h) and there is a saddle point 
at a symmetric linear configuration (-Doofi). The inclu- 
sion of nonadditivity deepens the linear and (especially) 
triangular trimer wells, but does not change the overall 
conclusions. 

The potential energy surfaces are such that barrier- 
less atom-exchange reactions can occur. All configura- 
tions of the trimer M3 are lower in energy than the sep- 
arated atom + diatom collision partners (and products). 
Once the collision complex M3 has been formed, any one 
of the three atoms can depart to form products. Even 
when the products are indistinguishable from the reac- 
tants, all three of these "arrangement channels" must 
be taken into account in the collision dynamics. This 
has two major consequences. First, a reactive scattering 
approach (rather than an inelastic scattering approach) 
must be used in the dynamics calculations. Secondly, 
the scattering calculations must be done on a fully three- 
dimensional potential energy surface. This is not the 
case for non-reactive scattering, where the products and 
reactants are confined to just one arrangement, and two- 
dimensional calculations (with the diatomic bond length 
kept fixed at the dimer equilibrium geometry) or quasi- 
three-dimensional calculations (where the diatomic bond 
length is varied only slightly around the dimer equilib- 
rium geometry) are usually used. 

The major focus of interest is in collisions that release 
kinetic energy and thus lead to trap loss. These are typ- 
ically vibrational relaxation processes of the type 

M 2 (v)+M — >M 2 (v'<v)+M. (4) 
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TABLE I: RCCSD(T) values of r c (A), r sp (A), V^in = -D e 
(cm -1 ), V sp (cm -1 ), and V3 (cm -1 ) for spin-polarized alkali 



dimers and trimers. 




Dimer 


Trimer D 


ih 


Trimer D 




r c 




r c 


Vini 11 


v 3 


' sp 


V Bp 


v 3 


Li 


4.169 


-334.046 


3.103 


-4022 


-5260 


3.78 


-968 


-354 


Na 


5.214 


-174.025 


4.428 


-837 


-663 


5.10 


-381 


-27 


K 


5.786 


-252.567 


5.084 


-1274 


-831 


5.67 


-569 


-52 


Rb 


6.208 


-221.399 


5.596 


-995 


-513 


6.13 


-483 


-15 


Cs 


6.581 


-246.786 


5.992 


-1139 


-562 


6.52 


-536 


-32 



However, if the three atoms are identical it is not possi- 
ble to distinguish between inelastic processes and reactive 
processes. We therefore use the general term quenching 
to describe collisions that produce a change in the vi- 
brational (or rotational) quantum number and release 
kinetic energy. 



POTENTIAL ENERGY SURFACES 
Nonadditive forces 

Accurate quantum scattering calculations require ac- 
curate potential energy surfaces. Higgins et al. [olf 
showed in 2000 that the quartet state of Na3 exhibits 
strong non-additive forces that increase the well depth 
of the equilateral trimer by 59% and decrease the Na- 
Na bond length from 5.2 A in the triplet dimer to 4.4 
A in the quartet trimer. We therefore carried out a 
systematic study 01 to investigate such effects for the 
whole series of homonuclear alkali-metal trimers. Ab ini- 
tio electronic structure calculations were performed using 
a single-reference restricted open-shell variant [H^] of the 
coupled-cluster method 63] with single, double and non- 
iterative triple excitations [RCCSD(T)]. Medium/large- 
size basis sets were used for the alkali metal atoms as 
described in ref. I47I and the full counterpoise correction 
of Boys and Bernardi [|34| was employed to compensate 
for basis-set superposition errors. All the ab initio cal- 
culations were performed using the MOLPRO package 
[65|. 

The results are summarized in Tabled which shows the 
equilibrium bond lengths r e and potential depths V m - ln 
for alkali-metal dimers and equilateral (-D3/1) trimers, to- 
gether with the corresponding quantities r sp and V sp for 
the linear (-Dooh) saddle points. The three-body non- 
additive contributions V3 are also given. Fig. ^ shows 
the additive and nonadditive potential energy curves for 
D^h geometries. 

All the trimers show quite strong nonadditive effects. 
The quartet trimers all have equilibrium interatomic dis- 
tances (at Z?3/j geometries) that are substantially shorter 
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FIG. 1: RCCSD(T) interaction energies of spin-polarized al- 
kali trimers at D 3 h geometries (a) full potentials including 
non-additive contributions; (b) additive potentials. Repro- 
duced from Soldan et al. 14711 - 



than those of the triplet dimers, by an amount that de- 
creases down the series from 1.07 A in L13 to 0.59 A in 
CS3. The trimer potentials are all correspondingly deeper 
than pairwise sums of dimer potentials, by a factor of 1.3 
to 1.5 for the heavier alkali metals (Na to Cs) but a factor 
of more than 4 for Li. 

The size of the nonadditivity is at first sight quite sur- 
prising in chemical terms. It contrasts with the situation 
for the rare gas trimers, where the non-additive contri- 
butions are only 0.5% to 2.5% j6|| and produce a 
weakening of the binding at equilateral geometries rather 
than a strengthening as in Table [fl 

The interaction potentials can be decomposed into 
self-consistent field (SCF) and correlation contributions. 
For the triplet alkali dimers, as for the rare gas dimers, 
the SCF potentials are repulsive and the main attrac- 
tive forces arise from interatomic correlation (dispersion). 
However, this similarity does not extend to the trimers. 
For the rare gases, most of the nonadditivity comes from 
the dispersion interaction. The leading long-range term 
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in this is the Axihod-Teller-Muto (ATM) triple-dipole 
term |68l l69| , which is repulsive near equilateral config- 
urations but attractive near linear configurations. For 
alkali metal atoms, by contrast, there is a large attrac- 
tive contribution to the non- additive energy that exists 
even at the SCF level. This occurs because the alkali 
metals have vacant np orbitals that lie relatively close to 
the ns orbitals. 

The question then arises why the p orbitals contribute 
so strongly for the alkali metal trimers but not the 
dimers. Soldan et al. |47| carried out a natural orbital 
analysis for equilateral triangle geometries, and consid- 
ered the contribution from radial p orbitals (pointing to- 
wards the centre of the triangle) and from tangential p 
orbitals (pointing around the ring). The radial p or- 
bitals can form bonding and antibonding molecular or- 
bitals (MOs) of the same symmetry as those formed from 
the ns orbitals (a[ and e'), while the tangential p orbitals 
can form a' 2 and e' MOs. The sets of MOs of the same 
symmetry interact, lowering the energy of the occupied 
MOs and contributing to bonding. Soldan et al. showed 
that the dominant contribution to the trimer bonding is 
from the tangential p orbitals. In chemical terms, this 
is essentially sp hybridization. This mechanism does not 
occur for the alkali metal dimers, because in that case 
the "tangential" p orbitals form molecular orbitals of 7r 
symmetry that cannot mix with a orbitals. 

Global potential energy surfaces 

It is relatively straightforward to generate potential en- 
ergy surfaces for systems such as the quartet alkali metal 
trimers by carrying out ab initio electronic structure cal- 
culations at grids of points that sample the configura- 
tion space. We have carried out such calculations for Li 3 
[Hli3 and K 3 jH3, using RCCSD(T) calculations as de- 
scribed above for the linear and equilateral geometries. 
Colavecchia et al. [7(J have also carried out calculations 
for Li 3 , using a more complete treatment of the valence 
electron correlation (full configuration interaction) but 
without correlating the core electrons. 

A much more difficult problem is to generate a global 
potential energy surface from a set of points. For quan- 
tum dynamics calculations, it is very important to repre- 
sent the potential-energy function smoothly and without 
oscillations between ab initio points. If the resulting po- 
tential is to be capable of representing all the properties 
of experimental interest (including atom-atom scattering 
lengths, dimer and trimer bound states, atom-diatom col- 
lisions and 3-body recombination) , then it is very impor- 
tant that the potential should dissociate properly into all 
possible sets of products (atom + diatom and 3 separated 
atoms) with the correct long-range behaviour. 

There are several coordinate systems that can be used 
for triatomic systems, including hyperspherical coordi- 



nates, Jacobi coordinates, and bond-length coordinates. 
These are by no means equivalent for interpolation pur- 
poses. In particular, grids of points in hyperspherical co- 
ordinates tend to include points in which 2 atoms lie very 
close together, which hinders interpolation because poly- 
nomials with very high localised maxima tend to have 
oscillations in other regions. Jacobi coordinates suffer 
from the same problem, and also do not allow the full 
3-body exchange symmetry to be introduced in a natu- 
ral way. Fortunately, there is no need to represent the 
potential energy surface in the same coordinate system 
as is used in the dynamical calculations. We therefore 
chose to carry out electronic structure calculations on a 
grid of points in bond- length coordinates {r\2, ^23, 7*13). 
We use the shorthand (r) for this to simplify notation. 

In order to represent the atom-diatom dissociation lim- 
its correctly, it is essential to use a long-range representa- 
tion in which the triatomic potential is decomposed into 
a sum of additive and non-additive contributions, 

3 

V(r) = W(r«) + V 3 (r). (5) 

i<j 

Provided Va(r) — > when any two of the atom-atom 
distances become infinite, this guarantees that the cor- 
rect diatomic potential is recovered in the atom-diatom 
limit. However, low-energy scattering is very sensitive 
to long-range forces, so we also require that the atom- 
diatom dispersion coefficients and their anisotropics are 
correctly reproduced. This requires careful treatment of 
the long-range part of V3 (r) . 

An important point is that, for a pairwise- additive po- 
tential in which the dimer potentials have the correct 
long-range form — CsR~ 6 , the atom-diatom C§ coeffi- 
cient is isotropic (independent of Jacobi angle 9). The 
anisotropy of the atom-diatom C§ coefficient comes en- 
tirely from non-additive forces. Cvitas et al. have 
therefore investigated the relationships between 3-body 
dispersion coefficients in the atom-diatom and atom- 
atom-atom representations and derived formulae relat- 
ing the atom-diatom C§ and C$ coefficients and their 
anisotropics to three-body coefficients arising from triple- 
dipole |6a,[6^, quadruple-dipole 0j[z3 and higher-order 
multipole |73J,|7J] terms. 

A variety of representations can be used for interpolat- 
ing dimer potentials. We have used the reciprocal-power 
reproducing; kernel Hilbert space method (RP-RKHS) 
[Til Ffij | . With an appropriate choice of parameters , 
this gives a potential that has the correct Cq and Cs 
dispersion coefficients. 

At short range, different approaches were needed for 
K3 and for Li 3 . K 3 is representative of the heavier alkali 
metals (Na to Cs) in that the nonadditive potential is 
smaller than the potential itself for most configurations. 
It is thus convenient to use Eq. [S] directly with a global 
representation of V3 (r) . The details of the procedure are 



given in refs. l5ll and 152. but in brief: 

• a quantity ^'(r) is defined by subtracting out 
damped versions of the triple-dipole and dipole- 
dipole-quadrupole terms from Va(r), 



y 3 '(r)=y 3 (r)- V^ p (v) + V^(v) 



(6) 



• a function 
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(7) 

is denned to eliminate the quadruple-dipole contri- 
bution to V^r); 

• the function V^"(r) = g(r) x V 3 '(r) is interpolated 
using 3-dimension RP-RKHS interpolation. V 3 "(r) 
is suitable for this (but V3 (r) and V^(r) are not) 
because Vg"(r) takes a product form at long range, 
constant x ^£2 ^23 '31 i 

• V^r) is then rebuilt from V^(v) at each interpo- 
lated point, 



V 3 (v) 



5(r) 



+ 



Vfl^U+V™*,® ■ (8) 



The resulting potential energy surfaces for spin-polarized 
K3 jHJ are shown in Fig. [5] Note that the depth at D^h 
geometries is rather more than twice that at linear ge- 
ometries, whereas pairwise additivity would give a factor 
of 1.5. 

Li 3 required a different procedure, because in this case 
the potential minimum for the trimer occurs at a dis- 
tance that is high on the repulsive wall for the dimer. 
Because of this, Eq. [S] would represent the interaction 
potential in this region as a difference between two very 
large quantities. Nevertheless, at long range a decom- 
position according to Eq. © is essential. Under these 
circumstances, it is best to fit the ab initio points di- 
rectly to obtain a short-range function VsR.( r ) without 
imposing the correct long-range behaviour. A switching 
function 5(r) is then used to join this onto the correct 
long-range form, 



V(r) = S(r)V SR (r) + [1 - S(r)]V LR (r). 



(9) 



The switching function is 1 at short range but switches 
smoothly to zero at long range. The long-range form 
Vlr(i") is designed to be valid when any of the atom- 
atom distances is large, and the procedure used to build 
in the correct three-atom and atom-diatom dispersion 
coefficients is described in ref. l52l 

A further complication arises for Li 3 because in this 
case there is a second potential energy surface involved. 
As described above, the quartet state that correlates with 
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FIG. 2: Cuts through the K3 quartet surface in valence coor- 
dinates. Upper panel: cut for a bond angle of 60°, showing 
the global minimum at —1269 cm -1 and 5.09 A. Lower panel: 
cut at collinear geometries; the collinear minimum is at —565 
cm -1 and 5.68 A. Contours are labelled in cm" 1 . Reproduced 
from Quemener et al. |5l| 



3 ground-state ( 2 S) atoms has 4 A' symmetry ( 4 S + at lin- 
ear geometries). The second state correlates with Li( 2 S) 
+ Li( 2 S) + Li( 2 P) and has 4 II symmetry at linear ge- 
ometries. It can cross the 4 E + state at linear geome- 
tries because it has different symmetry. However, at 
non-linear geometries the 4 I1 state splits into 4 A' and 
4 A" components, and the two 4 A' states mix and can- 
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r, / A 

FIG. 3: Cuts through the L13 quartet surface in valence co- 
ordinates. Upper panel: cut for a bond angle of 170°, show- 
ing the double-minimum structure due to avoided crossings 
at near-linear geometries. Lower panel: cut at collinear ge- 
ometries, showing the seam of conical intersections (red line) 
between the 4 E + and 4 IT states. Reproduced from Cvitas et 
al. [23. 



not cross. There are thus conical intersections at linear 
geometries, and the avoided crossings produce double- 
minimum structures at nonlinear geometries as shown in 
the upper panel of Fig. The line of conical intersec- 
tions produced a seam in the potential surface as shown 
in the lower panel. 



Conical intersections exist for both symmetrical and 
unsymmetrical linear geometries in quartet Li3, but for 
strongly unsymmetrical geometries the conical intersec- 
tions are high on the repulsive wall and will not af- 
fect the dynamics. However, for near-symmetric linear 
geometries the seam dips to an energy of about —100 
cnr 1 (relative to the energy of three free atoms) at 
t\i = r23 ~ 3.1 A. It is thus close to the inner turning 
point for low-energy collisions between Li and Li2, and 
may have significant consequences for the chemical dy- 
namics. The conical intersection has subsequently been 
characterized in more detail by Brue et al. [78| • 

QUANTUM DYNAMICS CALCULATIONS 
Methodology 

As described above, alkali metal atom + diatom colli- 
sions require a fully reactive scattering treatment, and at 
the energies of interest for cold molecules it is essential 
to handle relative translation as well as internal motions 
quantum-mechanically. Quantum reactive scattering cal- 
culations [Tij can in general be performed using either 
time- dependent or time-independent treatments. In re- 
cent work at higher energies, time-dependent treatments 
based on wavepacket dynamics have been becoming in- 
creasingly popular [sjj- At ultralow energies, however, 
the time evolution of a wavepacket is very slow, and prop- 
agating it until it reaches the asymptotic region requires 
an impractical number of time steps. Furthermore, it is 
difficult to converge wavepacket calculations at very low 
scattering energies, because the wave packet needs to be 
very broad. We therefore chose to describe ultralow en- 
ergy scattering in a time-independent formalism. 

There are many variants of time-independent reactive 
scattering theory. However, some of them are unsuit- 
able for the alkali metal trimers. As described above, the 
strong non-additive three-body interactions for the al- 
kali metal trimers make the atom-atom distances at the 
trimer equilibrium geometry much shorter than those for 
the dimers. The vibrational wavefunctions of a free dimer 
are centred around its equilibrium bond length and are 
very small at the distances that correspond to the trimer 
equilibrium. Because of this, the free dimer wavefunc- 
tions do not form a good basis set for expanding the 
scattering wavefunctions in the region of the trimer equi- 
librium, where the actual atom exchange takes place. 
This precludes the use of standard reactive scattering 
packages such as the ABC program |sjj , which uses such 
basis functions and has been widely used in studies of 
reactions such as F+H2 jUil 13 and 0( 3 P)+H 2 
at ultralow energies. 

We thus chose to use a scatteringformalism based on 
hyperspherical coordinates p,9,4> |86| : p is the hyperra- 
dius, which describes the size of the triangle formed by 
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the three atoms, while 9 and <fi are hyperangles, which 
describe the shape of the triangle. Hyperspherical ap- 
proaches do not use free-diatom functions as a basis set 
in the region of the trimer minimum. Instead, they define 
an adiabatic basis set by solving a fixed-p Schrodinger 
equation on a grid of values of the hyperradius p. There 
are several different hyperspherical approaches available. 
The approach developed for reactive scattering by Pack 
and Parker [&3> 0] solves the angular problem using a 
finite-element approach in adiabatically adjusting prin- 
cipal axis hyperspherical (APH) coordinates and then 
solves the resulting radial coupled equations by propa- 
gation. The approach developed by Esry and coworkers 
[Hit and used extensively for 3-body recombination 
in cold gases solves the angular problem in 

slightly modified Smith- Whitten coordinates using 
basis splines and then handles the radial problem with 
a finite clement approach. Both these methods use an 
optimised nonuniform grid in the hyperangles. However, 
we chose to use an alternative approach developed by 
Launay and LeDourneuf j^Dl > which has been applied ex- 
tensively to chemical reactions such as N( 2 D)+H2 |9l| 
and 0( 1 D)+H 2 [§2 at higher energies. 

In the approach of Launay and LeDorneuf, the config- 
uration space is divided into inner and outer regions, and 
the boundary between them is placed at a distance (hy- 
perradius) such that couplings due to the residual atom- 
diatom interaction can be neglected outside the bound- 
ary. This distance is typically p — 45 to 60 a for the al- 
kali metal systems. In the inner region, the wavefunction 
for nuclear motion is obtained by propagating a set of 
coupled differential equations using a diabatic-by-sector 
algorithm. The angular basis set is obtained by diag- 
onalizing a fixed-hyperradius reference Hamiltonian in 
a primitive basis set of pseudohyperspherical harmon- 
ics. In the outer region, the wavefunction is expanded 
in a basis set of diatom vibration-rotation functions ex- 
pressed in Jacobi coordinates |93|. The wavefunctions 
in the outer region are computed by inwards integration 
of regular and irregular solutions of an uncoupled radial 
Schrodinger equation which includes the isotropic part of 
the atom-molecule interaction. Matching between wave- 
functions for the inner and outer regions yields the scat- 
tering S-matrix. Elastic and inelastic cross sections are 
then calculated using standard formulae pSlj . 

The size of the basis set required for convergence de- 
pends strongly on the masses involved and the depth 
of the potential energy well. Before our calculations 
on alkali metal systems, nearly all quantum scattering 
calculations had been on systems containing only one 
or occasionally two non-hydrogen atoms. Calculations 
on alkali metal atoms require much bigger angular ba- 
sis sets, though fortunately at low energies only a few 
partial waves (values of J, the total angular momentum 
excluding spin) are needed. Even for Li + Li2 pjjj ]. which 
is relatively light, the number of adiabatic angular func- 



tions retained in the coupled equations ranged from 97 
for J = to 827 for J = 10, while for K + K 2 [Hl| the 
range was from 250 for J = to 1411 for J = 5. 

One problem with hyperspherical methods is that the 
diatom functions become more localised in hyperangular 
space as p increases. Because of this, the number of hy- 
perspherical harmonics needed for convergence increases 
with p: for K + K 2 it varied from 867 functions at small 
p to 6625 functions at the matching distance. It is the 
calculation to build the adiabatic angular basis that dom- 
inates the computer requirements (180 hours on an IBM 
Power4 P960 for K + K 2 ). 

A major advantage of the hyperspherical harmonics is 
that boson or fcrmion symmetry is very easy to impose. 
The complete nuclear permutation group for a system 
with three identical nuclei is 53. To satisfy the Pauli 
principle, the total wavefunction must have A\ symmetry 
for bosonic nuclei or A2 symmetry for fermionic nuclei. 
The total wavefunction is in general a sum of products 
of electronic, nuclear spin and nuclear motion parts. For 
three atoms in their spin-stretched states, the nuclear 
spin wavefunction is totally symmetric. For such states, 
collisions take place entirely on the quartet surface, for 
which the electronic wavefunction has A 2 symmetry. Bo- 
son or fermion symmetry can thus be imposed by select- 
ing pseudohyperspherical harmonics to give the correct 
symmetry for the wavefunction for nuclear motion. The 
adiabatic states in each sector are obtained by a varia- 
tional expansion on a basis of hyperspherical harmonics 
with A\ symmetry for bosonic atoms (with fermionic nu- 
clei) and Ai symmetry for fermionic atoms (with bosonic 
nuclei) . 

All the calculations were carried out on the quartet 
trimer surfaces, so are appropriate for collisions of spin- 
polarised atoms and molecules. However, the basis func- 
tions used for the quantum dynamics calculations did not 
explicitly include electron spin. In a more complete treat- 
ment, the rotational quantum number n for the triplet 
dimer would couple to its spin s = 1 to give a resultant j. 
When spin is neglected, however, there is no distinction 
between n and j. The splittings between levels of the 
same n but different j are in any case very small for the 
alkali metal dimers. 

Homonuclear molecules 

Quantum dynamics calculations have been carried out 
for the homonuclear collisions Li + Li 2 (t> = to 3) j3{J 

tNa + Na 2 (w = to 3) E3 and K + K 2 {v = 1) 
. For the Li and K systems, calculations were carried 
out for both bosonic and fermionic isotopes. 

The results for Li are typical. Figs.^]and[S]show elastic 
and inelastic cross sections for bosons ( 7 Li) and fermions 
( 6 Li) respectively. The elastic and inelastic s-wave cross 
sections for 7 Li are compared directly and extended to 
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FIG. 4: Elastic cross sections (upper panel) and inelastic 
cross sections (lower panel) for 7 Li + 7 Li2(vi = 1, j t — 0), 
with contributions from individual partial waves (total angu- 
lar momentum J, excluding spin). The vertical lines indicate 
centrifugal barrier heights for I > 1. Reproduced from Cvitas 
et al. |35|. 
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FIG. 5: Elastic cross sections (upper panel) and inelastic 
cross sections (lower panel) for 6 Li + 6 Li2(uj = 1, ji = 1), 
with contributions from individual partial waves (total angu- 
lar momentum J, excluding spin). The vertical lines indicate 
centrifugal barrier heights for I > 1. Reproduced from Cvitas 
et al. 1351 . 



lower energy in Fig. It may be seen that at very low 
energies (below 100 /iK for Li) the elastic cross sections 
become independent of energy whereas the inelastic cross 
sections are proportional to E-^JJ 2 . This is in accordance 
with the Wigncr threshold laws, which state that at very 
low energy the partial cross sections (contributions from 
a single partial wave I) for elastic and inelastic scattering 
vary as 

<4i - E% n ; a\ Bel ~ E [ - n ] /2 . (10) 

For a long-range potential proportional to R~ 6 , as for 
neutral atom-diatom scattering, there is an /-independent 
term that dominates the threshold law for high I so that 
°il ~ ^kin for / > 2 0. 

It may be seen from Figs. 0] and [S] that below E^ n = 



100 /iK the cross sections are completely dominated by 
the I = term, which corresponds to total angular mo- 
mentum J = for bosons but J = 1 for fermions (be- 
cause J = l+j and the lowest rotational level of a triplet 
fcrmion dimer is j = 1). The energy-dependence of the 
inelastic rate coefficient is fctTmoi where k = (2£ , km/A t ) 1/ ' 2 > 
so that the inelastic rate is independent of energy in this 
region. 

Above E'kin = 100 /iK, the I = contribution starts to 
deviate from the Wigner limit and higher partial waves 
start to contribute. The point at which this happens 
depends somewhat on mass, and is closer to 10 for K 

+ k 2 H2. 

As described in Section 2, it is known experimentally 
that fermion dimers produced by Feshbach resonance 
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FIG. 6: Elastic and inelastic s-wave (J = 0) cross sections 
for 7 Li + 7 Li2(vi = 1, ji = 0). The inset shows the real and 
imaginary parts of the complex scattering length. 



tuning are very much stabler than boson dimers when 
the scattering length is large. This stability was cru- 
cial in the production of molecular Bose-Einstein con- 
densates of 6 Li 2 H,@ and 40 K 2 [Oj. Petrov et al. |H 
explained the stability in terms of the requirements of 
Fermi-Dirac statistics. However, their derivation is valid 
only for long-range molecular states. A very important 
question is whether the stability persists for low-lying vi- 
brational states of fermion dimers. 

Cvitas et al. (3^ carried out quantum dynamics calcu- 
lations for spin-polarized collisions of ultracold homonu- 
clear Li + Li 2 collisions for both the bosonic ( 7 Li) and 
fermionic ( 6 Li) cases. The results shown in Figs. 0] and 
correspond to vibrational quenching rates for v = 1 of 
fc incl = 5.6 x 10~ 10 cm 3 s" 1 for bosons and 2.8 x 10~ 10 
cm 3 s _1 for fermions. However, the apparent difference 
of a factor of 2 is not significant: there were differences 
of up to a factor of 8 between inelastic cross sections for 
different initial values of v and j. Cvitas et al. found no 
systematic difference between inelastic rates in the boson 
and fermion cases, even when the atom-atom scattering 
length was adjusted to be large and positive. This has 
important consequences for efforts to transfer the dimcr 
population from Feshbach resonance states to the vibra- 
tional ground state, v — 0: it will be necessary to carry 
out the process either in a single step or quickly enough 
that the molecules do not spend enough time in interme- 
diate states to undergo inelastic collisions. 

The effective potential for a partial wave with I > is 
governed at long range by the centrifugal and dispersion 
terms, 



V\R) 



h 2 i(i + i) 

2/iR 2 



(11) 



where Cq is the atom-diatom dispersion coefficient. 
There is thus a centrifugal barrier at a distance 



6/iC 6 



h%l + l) 



1/4 



with height 



VI = 



h 2 i{i + i) 



3/2 



(54C 6 ) 



-1/2 



(12) 



(13) 



The resulting barrier heights are included in Figs. 0] 
and [5J The first vertical line corresponds to the I = 1 
partial wave and so on up to the I = 7 partial wave. It 
may be seen that each partial cross section has a maxi- 
mum at an energy slightly higher than the corresponding 
V^' lax . At collision energies below the centrifugal barrier, 
the partial cross sections for each I follow Wigner laws 
given by Eq. I|10|) . Above the centrifugal barrier, the in- 
elastic probabilities come close to their maximum value 
of 1 and the cross sections vary as E~ x because of the 
~kr 2 factor in the expression for the cross section. 



Capture model outside the ultracold regime 

At high collision energy, when several partial waves are 
involved, the total inelastic rate coefficient can be com- 
pared with that given by the classical Langevin capture 
model (93, based on the idea that every collision that 
crosses the centrifugal barrier produces inelasticity. This 
gives 
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This rate coefficient is shown as a function of collision 
energy for Li + Li 2 in Fig. [7] and compared with the 
quenching rates for bosons and fermions initially in v = 1 
and 2 35] . It may be seen that the full quantum result 
approaches the Langevin value at collision energies above 
about 10 mK. Similar behaviour is seen for K + K 2 at 
collision energies above 0.1 mK jiSlj . 



Product rotational distributions 

The vibrational spacings of the alkali metal dimers are 
much larger than their rotational spacings, so that many 
rotational levels are energetically accessible in collisions 
that cause vibrational relaxation. For Na 2 , for example, 
rotational levels up to j = 20 are energetically accessi- 
ble at the energy of the v = 1 state (23.5 cm -1 ). All 
accessible levels are populated in the products, subject 
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FIG. 7: Total inelastic rate coefficients for collisions of Li 
with Li 2 (v = 1 and 2, with j = for bosons and j — 1 for 
fermions). Reproduced from Cvitas et al. |3rI ]. 
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FIG. 8: Rotational distributions for 23 Na + 23 Na 2 (v = 1) 
at collision energy 10 -4 K. The label j' is the final rotational 
quantum number of 23 Na2 (V = 0). Reproduced from Soldan 
et al. El. 



to symmetry restrictions (only even-j levels for bosons 
and odd-j levels for fermions). The product rotational 
distributions for Na + Na2 (v = 1) at 10 -4 K are shown 
in Fig. [S] There are three clear maxima in the distribu- 
tion, at j = 4, 12 and 18. The oscillations probably arise 
from a rotational rainbow effect ^96J - The structure 
is similar to that observed in vibrational predissociation 
of Van der Waals complexes 98] . In a classical impulsive 
model, the energy released from Na2 vibration is partly 
retained in relative translation and partly converted into 
Na2 rotation. The angular momentum imparted to the 
Na2 molecule is zero if the energy is released at a linear 
or T-shaped geometry, but large around 9 — 45°. In this 
model, the oscillations arise from interference between 
classical trajectories on either side of the maximum. 

The rotational distributions become constant at low 
energies but show structure above 10~ 4 K. This is shown 
for K + K 2 (v = 1) in Fig. El In this case levels up to 
j — 24 are energetically accessible. Once again there is an 
oscillatory structure in the product state distributions. 




FIG. 9: Rotational distributions for 39 K + 39 K 2 (v = 1) as 
a function of the collision energy. The label j' is the final 
rotational quantum number of 39 K 2 (V = 0). Reproduced 
from Quemener et al. |5l|. 



Potential sensitivity 

The sensitivity of the cross sections to details of the 
potential energy surface is of great importance. Soldan 
et al. |4fij showed that including the nonadditive part of 
the interaction potential changed both elastic and inelas- 
tic cross sections for Na + Na2 (v = 1) by more than 
a factor of 10. Quemener et al. 48J investigated this in 
more detail: they introduced a scaling factor A to multi- 
ply the nonadditive term, so that A = corresponds to 
a pairwise-additive potential and A = 1 to the full non- 
additive potential. They then investigated cross sections 



as a function of A for initial v = 1,2 and 3. They found 
that elastic and inelastic cross sections varied by a factor 
of 10 for v = 1 for variations of A as small as 0.01 either 
side of A = 1. However, the variations became consid- 
erably smaller for v = 2 and 3. Cvitas et al. have 
investigated similar effects for Li + Li2 at rather higher- 
energy resolution. For v = 0, the elastic cross sections 
show very sharp structure as a function of A as shown 
in Fig. E| caused by poles in the scattering length each 
time there is a bound state at zero energy. However, the 
structure for v > is weaker for both elastic and inelas- 
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FIG. 10: Dependence of the elastic cross sections for 7 Li + 
r Li2(fi = 0,ji — 0) at E = 0.928 nK on the scaling factor 
A of the nonadditive part of the potential. Reproduced from 
Cvitas et al. E3- 



tic cross sections as shown in Fig. Ill| This is due to a 
general effect discussed by Hutson in the presence 
of inelastic scattering, the poles in scattering lengths at 
the positions of zero-energy resonances are suppressed, 
and the suppression increases with the degree of inelastic 
scattering. 



Differential cross sections 

At very low energies, cross sections are completely 
dominated by the I = partial wave. Under these cir- 
cumstances, the scattering is completely isotropic and 
the differential cross sections are featureless. However, 
as the energy increases and higher partial waves start to 
contribute, angular structure appears. Low-energy scat- 
tering thus offers the opportunity to study the onset of 
angular behaviour in reactive cross sections. The way 
that the angular behaviour develops is shown for K + 
K 2 (v = 1) in Fig. ^] At 1 /zK the scattering is com- 
pletely isotropic, but for 100 /xK some angular structure 
arising from interference between I = and 1 is evident. 
At 0.1 mK partial waves up to I = 5 contribute and sev- 
eral peaks emerge. 



Heteronuclear molecules 

Heteronuclear molecules are particularly interesting 
because they offer the possibility of studying reactive 
collisions separately from inelastic collisions. Homonu- 
clear molecules that are formed in their lowest vibration- 
rotation state are stable to collisions. For heteronuclear 
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FIG. 11: Dependence of the total elastic (upper panel) and 
inelastic (lower panel) cross sections for 7 Li + 7 Li2(uj, ji = 0) 
for Vi — 1, 2 and 3 and E — 0.928 nK on the scaling factor 
A of the nonadditive part of the potential. Reproduced from 
Cvitas et al. E3. 



molecules, however, the situation is more complicated. 
Even molecules in their ground rovibrational states may 
not be stable against collisions. For example, the spin- 
polarized reaction 

6 Li 7 Li(v = 0, j = 0) + 7 Li -> 6 Li + 7 U 2 (v = 0, j = 0) 

(15) 

is exothermic by 1.822 K because of the difference be- 
tween the zero-point energies of the two dimers. How- 
ever, the process 



6t ;7 



Li'Li(w = 0, j = 0) + b Li 



'Li 



3 Li 2 (« = 0,j = l) 



(16) 

cannot occur at collision energies below 2.643 K because 
of the combined effects of zero-point energy and the need 
to form 6 Li2 in j = 1 or higher to satisfy fermion sym- 
metry requirements. 
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FIG. 12: Differential cross section for inelastic scattering of 
K + K 2 (v = 1) at 1 fjK (dotted line), 100 (dashed line) 
and 0.1 mK (solid line). Reproduced from Quemener et al. 
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The results obtained in Ref. 50 have important impli- 
cations. There is interest in producing a quantum gas of 
6 Li 7 Li in its ground rovibronic state in an ultracold mix- 
ture of 6 Li and 7 Li atoms. In order to stabilize the molec- 
ular cloud against two-body trap losses induced by the re- 
active process (|15f) . the remaining atomic 7 Li would have 
to be removed quickly after ground-state molecule pro- 
duction, so that just the two-species fermionic mixture 
of 6 Li 7 Li(v = 0, j = 0) molecules and 6 Li atoms is left in 
the trap. The 6 Li cloud could be removed as well, but 
it might be advantageous to keep it in the trap. Elastic 
s-wave collisions between fermionic 6 Li 7 Li molecules will 
be strongly suppressed, but low-energy collisions with 6 Li 
can result only in elastic scattering and might be used to 
achieve sympathetic cooling of the molecules. 



Further extensions 
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FIG. 13: Elastic and reactive s-wave cross sections for 7 Li + 
6 Li 7 Li(u = 0, j = 0). The inset shows the higher-energy cross 
sections on a log-linear scale, with axes in the same units as 
the main plot. Reproduced from Cvitas et al. |H(J. 



Cvitas et al. [50j have investigated the process ljT5|l and 
the resulting elastic and reactive cross-sections are shown 
in Fig. It may be seen that the reactive scattering 
dominates over elastic scattering below 10 fiK . The low- 
temperature reactive rate coefficient is only 4.7 x 10~ 12 
cm 3 s , which is reduced by a factor of about 50 from 
those typical for vibrational relaxation in the homonu- 
clear case. Cvitas et al. attributed the difference to the 
fact that there is only a single open channel for reaction 

03. 



Our work on the quantum dynamics of collisions of 
alkali metal dimers has so far been restricted in several 
ways. We have focussed on collisions of molecules in low 
vibrational states, for systems involving three chemically 
equivalent atoms. We have restricted ourselves to col- 
lisions of spin-stretched atoms and molecules, for which 
doublet electronic states do not contribute. We have ne- 
glected hyperfine structure, and worked in zero applied 
field. 

Extending the calculations to handle heteronuclear 
systems is relatively straightforward, though a consid- 
erable amount of work is needed to develop potential 
energy surfaces for each system of interest. Dynamical 
calculations are more expensive for heavier atoms and 
for systems of lower symmetry (and the calculations de- 
scribed here already push the limits of current comput- 
ers). Extending the calculations to higher vibrational 
states is also mostly a matter of computer time, though 
true long-range states very near dissociation may be diffi- 
cult or impossible to converge with our current scattering 
methods. 

Including the effects of nuclear spin and magnetic fields 
is a very difficult task, though an important one if we are 
to explore atom-molecule Feshbach resonances and use 
them to control molecular interactions in the same way 
as atomic interactions. Collisions of atoms and molecules 
in non-spin-stretched states will be particularly challeng- 
ing, because they will involve doublet surfaces as well as 
quartet surfaces, and for alkali metal trimers the dou- 
blet surfaces exhibit co nical intersections and geometric 
phase effects |l00t Il0l| that considerably complicate the 
dynamics. 
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CONCLUSIONS 

This article focussed on theoretical studies of collisions 
between spin-polarized alkali metal dimers and atoms, 
which are crucial in experiments designed to form ultra- 
cold molecules in low-lying vibrational states. Colliding 
dimers can undergo very fast barrierless chemical reac- 
tions. As a result, vibrationally excited molecules un- 
dergo very fast vibrational relaxation, with rates usu- 
ally in excess of fci ne i = 1CP 10 cm 3 s _1 . At temper- 
atures above about 1 mK, where several partial waves 
contribute, the rates are approximately given by a sta- 
tistical Langevin capture model. At lower temperatures, 
however, the reactions enter a regime governed by Wigncr 
threshold laws and a full quantum-dynamical treatment 
is essential to calculate the rates. In this regime the re- 
sults are very sensitive to details of the triatomic po- 
tential energy surfaces, though the sensitivity decreases 
for excited vibrational levels. Isotopically heteronuclear 
molecules can often undergo exothermic reactions even 
from their ground vibrational states, because of the dif- 
ference in zero-point energy between reactants and prod- 
ucts. 

Prospects for the future include the production of 
quantum-degenerate gases of ground-state molecules, 
which will be stable to collisions and offer a wealth of 
new possibilities for quantum control. Heteronuclear 
molecules are particularly interesting, because they can 
have substantial dipole moments in long-range states. 
Dipolar quantum gases offer a new range of novel prop- 
erties, and ultracold polar molecules also have potential 
applications in quantum computing and in studying fun- 
damental physical properties such as parity violation and 
the electron dipole moment. 
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